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ABSTRACT 



Context. The prediction of stellar angular diameters from broadband photometry plays an important role for different applications. In 
particular, long-baseline interferometry, gravitational microlensing, extrasolar planet transits, and many other observing techniques 
require accurate predictions of the angular size of stars. These predictions are based on the surface brightness-color (SBC) relations. 
Aims. Our goal is to calibrate general-purpose SBC relations using visible colors, the most commonly available data for most stars. 
Methods. We compiled the existing long-baseline interferometric observations of nearby dwarf and subgiant stars and the correspond- 
ing broadband photometry in the Johnson BV and Cousins R C I C bands. We then adjusted polynomial SBC models to these data. 
Results. Due to the presence of spectral features that depend on the effective temperature, the SBC relations are usually not linear 
for visible colors. We present polynomial fits that can be employed with BVR C I C based colors to predict the limb-darkened angular 
diameters (i.e. photospheric) of dwarf and subgiant stars with a typical accuracy of 5%. 

Conclusions. The derived polynomial relations provide a satisfactory approximation to the observed surface brightness of dwarfs and 
subgiants. For distant stars, the interstellar reddening should be taken into account. 

Key words. Stars: fundamental parameters, Techniques: interferometric 



1. Introduction 

The surface brightness-color (hereafter SBC) relations link the 
emerging flux per solid angle unit of a light-emitting body to its 
color. They express the Stefan-Boltzmann relation between to- 
tal brightness and effective temperature in measurable quantities 
such as angular diameter, magnitude and color. It basically as- 
sumes that temperature and bolometric correction can be trans- 
lated into a color measurement. These relations have many as- 
trophysical applications, such as the estimation of Cepheid dis- 
tances (through the Baade-Wesselink method), extrasolar planet 
transit studies or the characterization of microlensing sources. 
Different colors produce tighter or more dispersed relations, sim- 
ply because color also depends on other variables such as gravity 
or metallicity, at a different level depending on the adopted pho- 
tometric bands. It has been shown that a combination of visible 
and near-IR bands, such as (V-K), is probably optimal (see, e.g., 
Fouque & Gieren [1997l in the case of Cepheids). However, near- 
infrared photometry is not always available, for instance due to 
catalogue incompleteness or field crowding of IR atlases such as 
2MASS in the Galactic disk. As a consequence, SBC relations 
in the visible domain remain very useful. 

Although they are relatively similar for stars of different 
luminosity classes, the literature gives appropriate relations 
for Cepheids (Kervella et al. |2004cl >, giants (van Belle I1999I 
Nordgren et al. I2002I) . M giants and supergiants (van Belle 
et al. 1 1 9991 Groenewegen 2004), and dwarf or subgiant stars 
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(Kervella et al. 2004b, hereafter K04). Here we limit our discus- 
sion to dwarf stars and subgiants, and revise and extend K04's 
work, which was based on Johnson photometric bands and was 
limited to linear SBC relations. In the present Note, we cali- 
brate the SBC relations based on Cousins R c and I c photometric 
bands, of more common use than Johnson's R and /, and we ad- 
just non-linear SBC relations that are a better fit to the observa- 
tions than linear laws. In this process, we use exclusively direct 
interferometric angular diameter measurements, thus avoiding 
cross-correlations with previously calibrated SBC relations. 

2. Interferometric and photometric data 

We collected from the literature all the available interferomet- 
ric angular diameter measurements of dwarf and subgiant stars. 
From this list, we removed the stars for which the angular di- 
ameter was uncertain by more than 5%. The major sets of new 
interferometric measurements of dwarf and subgiant stars since 
K04 are from Berger et al. d2006l l , Bain es et al. (|5o08), Boyajian 
et al. (12008b and Kervella et al. ((2008), all four based on mea- 
surements obtained with the CHARA array (ten Brummelaar 
et al. 120051 1. The only star from Berger et al. ( 120061 ) that has 
an angular diameter more precise than 5%, GJ 15A, was re- 
moved from our list as it is a flare star, and its photometry 
is therefore uncertain. This was also the case for several other 
M dwarfs. Most of the candidate stars present variability at a 
certain level. In order to reject those stars for which the vari- 
ability amplitude could be large, we checked their variability 
status i n the GCPD (Mermilliod et al. [19971 ) and the GCVS 
(Samus 120081 . All stars listed in "yCas" or "UVCeti" classes 
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have been rejected. Low amplitude "5 Scuti", such as a Lyr or 
B Leo, "BYDra" like e Eri, GJ 570 A, GJ 699, 61 Cyg A, HR 
511, or possible "RS Cvn", such as 6 Eri, have been kept. As a 
remark, DY Eri (a "UV Ceti" type flare star) is not GJ 166 A, as 
stated in the SIMBAD database, but GJ 166 C, therefore we kept 
the former in our list. We also removed the known very fast ro- 
tating stars (v sin i > 100 km.s -1 ), including the pole-on rotator a 
Lyr (Aufdenberg et al. 2006; Peterson et al. 2006 ). The effective 
temperature of fast rotating stars is inhomogeneous over their 
apparent disk, they can be surrounded by circumstellar material 
(Kervella & Domiciano de Souza 2006 ) and their photospheres 
can be spectacularly distorted (Monnier et al. [2007). All these 
phenomena combine to bias their surface brightness compared 
to normal stars. Overall, our selection procedure resulted in the 
removal of about one third of the available interferometric mea- 
surements (Table [TJ. 

The limb darkened (LD) angular diameters were available 
in all cases from the original publication, except the measure- 
mentd of Sirius by Hanbury Brown et al. (119741 ) and Davis 
& Tango (1986), for which we applied the LD correction by 
Claret (120001 ). as discussed in Sec t. 5.4 of Kervella et al. d2004at . 
For Mozurkewich et al.'s ( 120031 ) measurements we considered 
the combined LD values listed by these authors. In all cases, we 
chose to keep the original LD values instead of correcting the 
uniform disk values by recomputing the limb darkening correc- 
tions, because there is a remarkable concensus in the literature 
to use the predictions from Kurucz's atmosphere models, ap- 
proximated using Claret's laws (Claret et al. 119951 Claret[2000i 
Claret 2000) to estimate the brightness distribution over the 
stellar disk. In fact, all the listed angular diameters were com- 
puted using these LD models, except the recent measurement of 
aCen B published by Bigot et al. (2006 ) who employed 3D hy- 
drodynamical models. Even in this case, the derived angular di- 
ameter was less than lcr away from the ID Kurucz model result. 
Claret (2008) recently showed that theoretical atmosphere mod- 
els do not perfectly reproduce the observed stellar intensity pro- 
files from transiting exoplanets and eclipsing binaries. However, 
these discrepancies are negligible for our purpose, considering 
the star-to-star residual dispersion with respect to the adjusted 
empirical laws (Sect. [3). 

The starting point for the photometry listed in Table [2]is the 
Hipparcos Input Catalogue (1993 ) as retrieved from SIMBAD, 
which lists V and B-V, and the Hipparcos and Tycho Catalogues 
( 1 19971 ), which list V-I. Notes in these catalogues give the origin 
of the listed values. When original photometry in Cousins sys- 
tem exists for V - R c and V - I c , it was checked from the on-line 
version of the Lausanne General Catalogue of Photometric Data 
d 19971 ), averaged over original sources according to their number 
of measurements, and finally combined to the V band Hipparcos 
photometry to get R c and I c magnitudes. More details are given 
in the Notes of Table HI 

3. Surface brightness-color relations 

Contrary to visible-infrared colors from K04, the dependence 
of the zero-magnitude limb-darkened angular diameter (ZMLD, 
defined for = 0), as a function of the color is non linear in 
the visible. For this reason, we selected polynomial SBC relation 
models of the form: 

log(#LD) = a + a\ C + a 2 C 2 + a 3 C 3 + a 4 C 4 - 0.2 m A (1) 

where C is the color of the star from two dereddened photometric 
bands among BVR C I C (e.g. B — V, V - R c ,---), the dereddened 



Table 1. Interferometric data. The error bars are listed in small 
characters as exponents. 
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magnitude of the star in one of the bands, and #ld the limb- 
darkened angular diameter, measured in milliarcseconds (mas). 
Note that this expression is independent of the distance of the 
star. For the empirical fit to the observations, we selected the 
smallest polynome degree that gave a satisfactory fit of all the 
data points. This led us to select a higher degree (four) for the 
B-V based polynomial fits. 

We would like to stress that the computed relations are valid 
only for a given range of colors, with only little margin for ex- 
trapolation. In particular, the SBC laws undergo steep variations 
for red objects. The results of the fits are presented in Fig.[T]for 
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Table 2. Photometric data. When the uncertainty on the photo- 
metric measurements was not available, we adopted an arbitrary 
+0.05 value. 
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B and V and in Fig.|2]for R c and I c . The polynomial coefficients 
and the domain of validity of each relation, are listed in Table [3] 
The residual dispersions ct(0 ld )/6'ld give the relative accuracy 
of the angular diameter predictions for each relation. 

4. Discussion 

4. 1 . Reddening corrections and metallicity 

When a SBC relation is applied to a star suffering some amount 
of absorption, its colors must be corrected for extinction before 
applying the relation. The amount of extinction is generally mea- 
sured by the color excess E(B - V), and the conversion to other 
photometric bands depends on the intrinsic color of the star and 
on the amount of extinction (Dean et al. I19781 ). For blue stars, 
the color excess ratios are typically E(V - R c ) = 0.60 E(B - V) 
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Fig. 1. Adjusted polynomial relations giving the zero-magnitude 
limb-darkened disk angular diameter (ZMLD 1) in the B and V 
bands as a function of different color combinations. The corre- 
sponding polynomial coefficients are listed in Table [3] 



and E(V —I c ) = 1.25 E(B-V). For red stars, the coefficients may 
increase by up to 10%. 

All the stars in our list are very nearby. In average, we do 
not expect a significant deviation of their metallicity compared 
to that of the Sun. The exoplanet stars observed by Baines et 
al. ([2008) probably present a slight overmetallicity, but the ef- 
fect on the dispersion of the adjusted SBC relations is taken into 
account in the dispersions stated in Table [3] As discussed by 
K04, the effect of metallicity on the visible-infrared SBC rela- 
tions is undetectable. For the visible SBC relations presented 
here, different metallicities will translate into more significant 
photometric color differences. However, we expect them to be 
small compared to the measured intrinsic dispersion of the rela- 
tions, at least over a +0.5 dex range around solar metallicity. 

4.2. Application to transits and microlensing 

The existence of exoplanetary transits in front of a star allows 
to retrieve its linear radius, as well as the planetary radius. 
However, when the parallax is not known accurately and/or the 
star is too faint for high accuracy angular velocity measurements 
(from high resolution spectroscopy), they are based on an a pri- 
ori radius for the star. The relations established in the present 
Research Note are useful in this context, as they allow to retrieve 
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Table 3. Polynomial coefficients of the adjusted surface brightness-color relations: log(#LD) = ao+#i C+CI2 C 2 +«3 C 3 +«4 C 4 — 0.2 
(Old expressed in milliarcseconds). N is the number of measurements used for the calibration. The listed C m j n and C max values 
correspond to the minimum and maximum color index of the stars included in the polynomial fitting procedure. They represent the 
limits outside of which the adjusted relations are unreliable and should not be employed. The standard deviation of the residuals of 
the fit in log scale is listed in the cr(log # LD ) column, and the corresponding relative uncertainties on the predicted angular diameter 
is given in the rightmost column, in percentage. 
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Fig. 2. Adjusted polynomial relations in the R c and I c bands. 



the angular size of faint stars for which only broadband photom- 
etry is available. As a test, we can apply the (7, V - I c ) relation 
to the nearby transiting exoplanet stars HD 189733, for which 
7 = 7.686 + 0.007, and (7 - I c ) = 0.93 + 0.01 (from Hipparcos). 



Due to the proximity of this star, we neglect the reddening. For 
these bands, the surface brightness relation is (TableO: 

log(6» LD ) = 0.4992 + 0.6895 (7-/ c ) - 0.0657 (V- I c ) 2 - 0.2 7(2) 

and we obtain 6» LD (HD189733) = 0.352 + 0.020mas. The at- 
tached 6% uncertainty contains the photometric uncertainties 
and the intrinsic dispersion of the [V, V - I c ) relation. Together 
with the Hipparcos parallax of this star n — 5 1 .94 + 0.87 mas, we 
obtain a linear photospheric radius of 0.728 + 0.044 R Q . This is 
in excellent agreement with the radius derived from HS T transit 
observations by Pont et al. ( |2"0"0"7j l of 0.755 + 0.01 1 R G . We can 
also compare the (V, V - I c ) angular diameter to the one com- 
puted from the infrared [V, V - K) linear relations calibrated 
by K04. The 2MASS catalogue (Skrutskie et al. 2006} lists 
K = 5.541 +0.021, giving LD (HD189733) = 0.372 + 0.007 mas. 
Although this value is within lcr of the (V, V - I c ) value, the un- 
certainty on this visible-infrared prediction is only one third of 
the visible version. A complete review of the properties of this 
system can be found in Torres, Winn & Holman (2008 ). 

Another important application of SBC relations is related 
to gravitational microlensing studies, when one needs to esti- 
mate the source radius from color-magnitude diagram (CMD) 
in the direction of the Galactic Bulge. Deredenned source mag- 
nitude and color are obtained from the relative position of the 
source in a CMD, with respect to the mean position of the Red 
Giant Clump (RGC). As the intrinsic position of the RGC is cal- 
ibrated from nearby stars by Hipparcos, the assumption that the 
source suffers the same amount of extinction as the clump avoids 
the need for an estimate of the extinction (see Yoo et al. 2004 
for details). Generally, the magnitudes are measured in the I c 
band, and the measured color is (V - I c ), because this is what 
OGLE measures. A caveat of K04's SBC relations is that they 
are attached to the Johnson R and / bands, and this has been 
a source of confusion in the past. Let's take the example of the 
recently discovered 3 Earth mass planet event MOA-2007-BLG- 
192 (Bennett et al. 2008 ). The source appears to be fainter than 
the RGC by 5.7 mag in I c , and redder by 0.07 mag in (V - I c ). At 
the source distance (7.5 kpc), the RGC is at 7 e0 (RGC) = 14.13 
and (7 - 7 C ) (RGC) = 1.04, so the source has 7 c0 (MOA 192) = 
19.84 + 0.24 and (7 - 7 C ) (MOA 192) = 1 . 1 1 ± 0.24 (see Bennett 
et al. 120081 for a detailed discussion of these values). Applying 
the appropriate (I c , V - I c ) SBC relation from Table [3] gives an 
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estimated angular diameter of # LD (MOA 192) = 0.99 + 0.14 yuas, 
which is resolved by the caustic crossing (or cusp approach). 

4.3. Comparison with other calibrations 

In their early publication, Barnes et al. (1 19781 1 obtained a lin- 
ear SBC relation F v = 3.964 - 0.333(5 - V) for 0.10 < 
(B-V) < 1.35, with a typical residual scatter of cr(F v ) = +0.025. 
Following K04, it can be expressed as log(6>LD) = 0.5134 + 
0.666 (B -V) - 0.2 V with a scatter of cr(log6> LD ) = +0.05, 
equivalent to <x(# LD ) = +12%. Over their common range of ap- 
plicability, the relative difference p in 6*ld with our polynomial 
relation (coefficients listed in Table [3} is 

p(B78) = [0 LD (B78) - LD (KO8)] /0 ld (KO8) = 3.6 + 3.5%. (3) 

We can also compare our (B, B - I) predictions (as an ex- 
ample) with the {V, V-K) relations calibrated recently by 
Beuermann (120061 hereafter B06). The Sy(V - K) relation of 
this author's Table 2 can be reformulated as: 

log6> LD =0.524 (V-K) 2 + 0.267 (V-K)- 1.221 10^-0.2 V(4) 

For the 34 stars in our sample, we obtain an average value of 
the relative difference with our calibration of p(B06) = 4.1 + 
5.2%. One should note that B06's calibration relies on several 
data sources, including indirect angular diameter estimates. The 
same comparison with K04's (V, V-K) relation gives p(K04) = 
0.7+4.6%. The present calibration of the visible SBC relations is 
therefore compatible with previously determined SBC relations 
using visible and visible-IR color indices. 

5. Conclusion 

We computed polynomial SBC relations in BVR C I C that can be 
used to predict the angular size of individual stars for which 
only broadband photometry is available. The visible-infrared lin- 
ear relations established by K04 usually provide more accurate 
angular size estimates. However, in many cases the JHK band 
magnitudes of field stars are not available (due to crowding in 
2MASS for instance), and the present BVR C I C relations will pro- 
vide a reliable photospheric angular diameter prediction with a 
typical uncertainty of m 5%. Their compatibility with visible-IR 
relations gives further confidence in their accuracy. 
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